Derivation of magnetic Coulomb's law for thin, semi-infinite solenoids 
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It is shown that the magnetic force between thin, semi-infinite solenoids obeys a Coulomb-type 
law, which corresponds to that for magnetic monopoles placed at the end points of each solenoid. 
We derive the magnetic Coulomb law from the basic principles of electromagnetism, namely from 
the Maxwell equations and the Lorentz force. 



I. INTRODUCTION 

A permanent magnet is an ensemble of microscopic 
magnetic moments which are oriented along the magne- 
tization direction. A magnetic moment can be modeled 
as a dipole, i.e., as a slightly displaced pair of magnetic 
monopoles with opposite polarities. It is analogous to the 
electric dipole. Another way of modeling is to consider 
each magnetic moment as a circulating current loop. In 
terms of far fields, the dipole model and the loop-current 
model give exactly the same magnetic field. The lat- 
ter model is more natural because there are no magnetic 
monopoles found so far and microscopic magnetic mo- 
ments are always associated with kinetic rotations such 
as orbital motions or spins of electrons. It also provides 
correct symmetries with respect to the time and space 
inversions. 

Normally we are interested in macroscopic quantities, 
which are obtained by coarse-graining of microscopic 
fields and source distributions^. When we coarse-grain 
the oriented ensemble of microscopic magnetic dipoles in 
a bar magnet, we have a magnetic north pole at one end 
and a south pole at the other end as shown in Fig. 1(a). 
Contributions of the magnetic charges inside are canceled 
out through the spatial average. 

On the other hand, when we coarse-grain the micro- 
scopic loop currents, we have macroscopic current which 
circulates around the bar along the surface as shown in 
Fig. 1(c). By mimicking the macroscopic current distri- 
bution by a coil, we have an electromagnet that is equiv- 
alent to the permanent magnet. In this model, the poles, 
or the ends of bar magnet play no special roles. Again 
the current model is much more reasonable because of 
the similarity with the equivalent electromagnet. 

The use of the notion of magnetic pole should be 
avoided as far as possible^ because of its absence in the 
framework of electromagnetic theory or in the Maxwell 
equations. Practically, however, magnetic poles are very 
convenient to describe the forces between permanent 
magnets or magnetized objects. The poles are consid- 
ered as an ensemble of monopoles and the forces be- 
tween poles are calculated with the magnetic Coulomb 
law, which is usually introduced just as an analog of the 
electric Coulomb law or as an empirical rule^. For logical 
consistency, we have to derive the magnetic Coulomb law 



from Maxwell's equations and the Lorentz force, none of 
which contain the notion of magnetic monopoles. 

The derivation of magnetic Coulomb's law was given 
more than forty years ago by Chen* and Nadeau^. In 
this paper a more detailed analysis based directly on the 
fundamental laws will be provided. The field singularity 
for infinitesimal loop currents, which plays crucial roles 
but was not mentioned in the previous works, will be 
treated rigorously. 



II. CURRENT DENSITY FOR A THIN 
SOLENOID 

For brevity, we introduce a scalar function Go and a 
vector function G\ of position r = (x, y, z): 
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We note that VG = -G x and V • G x = S 3 {r) hold, 
where <5 3 (r) = 6(x)5(y)5(z) is the three dimensional delta 
function. With these, the scalar potential for a point 
charge q placed at the origin is (f>(r) — (q/eo)Go(r), and 
the force acting on a charge q± at r\ from another charge 
q 2 at r 2 is F^ 2 = {qiqi/ £o)G 1 {r 1 - r 2 ). The Biot- 
Savard law can be expressed as dH — dC x Gi (r ) , where 



magnetic dipole 



current loop 
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(b) 



(c) 



FIG. 1: A permanent magnet is an ensemble of microscopic 
magnetic moments (b). If we consider each magnetic moment 
as a magnetic dipole, the corresponding macroscopic picture 
is two opposite magnetic poles at each end (a). If we adopt 
the current loop model for magnetic moment, the macroscopic 
picture consists of circulating current around the side wall (c). 
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dC = Idl is a current moment (current I times length 
dl) located at the origin. 

As shown in Fig. 2, a thin solenoid can be constructed 
as a stack of tiny loop currents at least in principle 6 . The 
current density for an infinitesimal loop current place at 
the origin is 

J m (r) = (-m x V)<5 3 (r) (2) 

where m is the magnetic moment for the loop current and 
its unit is (Am 2 ) (See Appendix). To form a solenoid we 
stack these tiny loop currents along a curve L a with a 
constant line density K a (loops per unit length). Each 
loop are aligned so that the direction of tangent of curve 
and that of m coincide. 

The current density distribution for a line segment 
dZ a ( || dm ) is 

dJ(r) = {-C a Ah x V)5 3 {r - r a ), (3) 

where C a = n & m ~ A m is the magnetic moment per unit 
length and characterizes the strength of the solenoid 7 , 
where m — m ■ dZ a /|dZ a | is the magnitude of m and r a 
represents the position of the segment dZ a . 

III. FIELD BY A THIN, SEMI-INFINITE 
SOLENOID 

Here we introduce an important formula. For a con- 
stant vector a. and a vector field V(r), we have 

(a x V) x V - a x (V x V) 

= (a • V)V - a(V • V). (4) 

With the Biot-Savard law and Eq. @, we can find 
the strength of magnetic field generated by a magnetic 
moment m placed at the origin as 

H m (r) = J dv'J m (r') xGx(r- r') 

= (-m x V) x G^r), (5) 



where dv' is a volume element at r' and V = d/dr'. 
Using Eq. Q, it can be rewritten as 

H m (r) = -(m- V)G 1 (r)+m6 3 (r), (6) 

where the relations V • G\(r) = 5 3 (r) and V x G\ = 
have been used. As shown in Appendix, Eq. ^ can also 
be derived from the Maxwell equations. 

It is well known that the field for an electric dipolc 
moment p is D(r) = — (p ■ V)Gi(r), which contains no 
delta-function terms unlike Eq. This means that in 
terms of near fields, the dipole and the current loop do 
not yield the same field. 

From Eq. the strength of the magnetic field created 
by a line segment dZ a is 

dfl(r) = (C a dZ a • V a )d(r - r a ) + C,dU 3 (r - r a ) 

(7) 

where V a = <9/<9r a . The integration along a curve L a 
yields the magnetic field created by the solenoid: 

H(r) = f dH = C a [Gx(r - n) - d(r - r 2 )] 

+ C a / ,5 3 (r - r a )dZ a , (8) 

where r 2 and T\ arc the start and the end points of 
L a . The second term of the right-hand side, which cor- 
responds to the magnetic flux confined in the solenoid, 
vanishes outside. For semi-infinite cases (r 2 = oo), we 
have 

H{r) = C & G l {r-r 1 )+C a , [ S 3 (r - r a )di a . (9) 

The first term of the right-hand side, which represents 
the magnetic field outside of the solenoid, is equivalent 
to the field for a monopole g a — /ioC a located at r\\ 

B{r)= jh G 1 {r-r 1 ) = i±^l s . (10) 
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The dimension of g a , 

9a = MoCa~-Am=^^Am = Vs = Wb (11) 
m m 

correctly corresponds to that for the magnetic charge. 
The magnetic flux g a confined along the solenoid fans out 
isotropically from the end point r\. As seen in Fig. 2, a 
thin, semi-infinite solenoid can be viewed as a magnetic 
monopole located at the end. 



IV. MAGNETIC FORCE ACTING ON A 
SEMI-INFINITE SOLENOID 

From Eq. @, we see that the Lorentz force acting on 
a tiny loop current m placed at r in a magnetic field B 
is 

F m = Jdv'J m (r') x B(r + r') 

= (m x V) x B(r). (12) 

Using Eq. (H|) and the conditions for magnetic field: V • 
B = (divergence-free) and V x (fj,Q 1 B) = (rotation- 
free), the expression can be modified as 



F m = (m- V)B(r), 



(13) 



which is suitable for line-integral. Equation (fT2"|) can also 
be modified as F m — V(m • B) with the divergence- free 
condition only. The rotation-free condition is satisfied 
only when the right hand side of the Maxwell- Ampere 
equation, V x H = J + dD/dt, vanishes. 

Thus the magnetic force acting on a line element dZb 
at »"b is 



dF = (C h dl h • V b )B(r b ), 



(14) 



where Cb = n^m and Vb = d/dr^. Integration along a 
curve Lb yields the total force acting on the solenoid; 



/ dF = C h I (dZ b • V b )B 
C h [B(r 3 ) - B(r 4 )] , 



(15) 



where r 4 and r 3 are the initial and end points of -Lb, 
respectively. For semi-infinite cases [B(r 4 = oo) = 0], 
we have 



F = C b B(r 3 ) = g h H(r 3 ), 



(16) 



with H — fJ.Q 1 B. This is equal to the force for a point 
magnetic charge t?b = MoCb placed at r 3 . 

It is surprising that the sum of the forces acting on each 
part of the solenoid can be represented in terms only of 
the magnetic field B(r 3 ) at the end point r 3 . This is 
because the partial force is proportional to the (vector) 
gradient of the magnetic field. 



V. COULOMB'S LAW BETWEEN TWO THIN, 
SEMI-INFINITE SOLENOIDS 



Now we can calculate the force (| L6[) on solenoid b by 
the field @ generated by solenoid a; 



C h B(r 3 ) = C h ^ H(r 3 ) 
MoCbC a Gi(r 3 - n). 



(17) 



This equation corresponds to Coulomb's law for magnetic 
charges g a at n and t/b at r 3 : 



F 3 ^i 



-Gi(r 3 - n) = -. (18) 



/'o 



4tt^ |r 3 - ri| 



These two equations are exactly the same but the for- 
mer is for the total (integrated) force between currents 
of solenoids and the latter is for the force between mag- 
netic point charges. 

It is interesting that the force is independent of the 
paths of either solenoids as long as their end points and 
strengths are kept constant. We should note that the 
derivation of the magnetic Coulomb law is not straight- 
forward. The properties of static, source-free fields must 
be utilized. 



VI. TORQUE ON A SEMI-INFINITE SOLENOID 

Now we are interested in the torque due to the inte- 
grated magnetic force on a thin, semi-infinite solenoid. 
The torque N m with respect to the origin exerted on a 
loop current m located at r is^ 

N m = J Av'{r + r') x [J m {r') x B(r + r')] 

F m + J dv' [B(r + r') x J m (r')} x r' 



= r x 
= + iV« 

- 1 Tfi. < * m. 



(19) 



where the first term corresponds to the torque due to the 
total force F m with respect to the origin and the second 
term corresponds to the torque with respect to the center 
of the loop current. Using Eq. (|T5|) . the former can be 
modified as 

=r x F m = r x (m ■ V)B 

= (m ■ V)(r x B) - \{m ■ V)r] x B 

= (m -V)(r x B) - m x B, (20) 

with (m ■ V)r = m. The latter can be simplified as 

N ™ = J dv 1 [B(r + r') x (-m x V)S 3 (r')] x r' 

= r'x[(mxV')x5(r + r')]| r , =0 

+ [B(r + r') x (m x V')] x r'| /=0 
= + [(B • V')m - (B • m)V] x r'\ r , =(j 
= mxB, (21) 
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where we have utilized V' x r' = with V' = d/dr' . . 
Finally, we have 

N m = N^+N^ = (m-V)[rxB(r)}. (22) 

The torque on a line element dZ b at r b is 

dN = C h (dl h - V b )[r b x B(r h )}, (23) 

and the integration along the semi-infinite solenoid L b 
terminated at r 3 yields 

N= [ dN = C h [ (dZ b • V b )[r b x B{r h )} 
J Lh Jh h 

= C h r 3 x B(r 8 ) = r 3 x F, (24) 

where F is the force (|16l) on the solenoid. Surprisingly, 
the torque exerted on a thin, semi-infinite solenoid coin- 
cides with that for a point magnetic charge gb = A<oC b 
place at r 3 . 

VII. DISCUSSION 

In this paper we only deals with thin solenoids. 
Solenoids with finite cross-section can be represented as a 
bundle of thin solenoids. Similarly a permanent magnet 
can naturally be modeled as a bundle of thin solenoids. 
The degree of magnetization of magnets or magnetized 
objects is characterized by the quantity called macro- 
scopic magnetization, M ~ A/m, which is defined as 
volume density of magnetic moments, m ~ A m 2 . The 
magnetization can also be represented as area-density of 
thin solenoids. We remember that a thin solenoid is char- 
acterized by the current moment, C~ Am. Bundling v 
solenoids per unit area (~ /m 2 ), we can create the mag- 
netization M = vC. 

In conclusion, we have shown that the force between 
thin, semi-infinite solenoids obeys the Coulomb law that 
is for the equivalent magnetic point charges placed at 
each end of solenoid. We have also shown that the torque 
exerted on the solenoid in a magnetic field coincides with 
that for a corresponding magnetic point charge. It is 
convenient to introduce magnetic charges or poles be- 
cause the magnetic Coulomb law can easily be applied 
for forces between them. But we should remember that 
without the justification given in this paper it is just a 
rough and ready method. 
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APPENDIX: CURRENT DENSITY AND FIELD 
OF A TINY CURRENT LOOP 

The current density @ for a tiny current loop can be 
derived as follows. We consider a parallelogram defined 



by a pair of small vectors, a and b. The center is located 
at the origin and current / is circulating along the edge. 
The current distribution on the four segments can be 
approximately represented as 

Ji(axb)(r) =1 [aS 3 (-b/2) + bS 3 (a/2) 

-aS 3 {b/2)-b5 3 (-a/2)] (A.l) 

In the limit of \a\, \b\ — > with m = I(a x b) being kept 
constant, it approaches 

J m (r) = I [-a(b ■ V)5 3 (r) + b(a ■ V)5 3 {r)] 
= [-I{a xb)x V}S 3 (r) 
= (-m x V)8 3 (r), (A.2) 

and Eq. ([2|) is obtained. The derivative of the delta func- 
tion V<5 3 (r) serves as a differential operator when it is 
integrated together with other functions as in ID cases: 
/ dxf(x)(d/dx)6(x) = — (d/ '/dx)(0) . For example, the 
magnetic force on m can be calculated as follows; 

F m = J dv(-m x V)(5 3 (r) x B(r) 

= J dv [-V5 3 (m ■ B) + m(V8 3 ■ B)} 

= [V(m • B) — m(V • B)] r=a 

= [(m x V) x B] (0). (A.3) 

Next we derive Eq. (|6|) from Maxwell's equations. The 
Ampere law for an infinitesimal loop current, V x H = 
J m can be written as 

V x (H(r) - md 3 (r)) = 0. (A.4) 

We know that a rotation-free field can be represented as 
a gradient of some scalar field </> m (r) as 

H{r) - md 3 (r) = -V^ m (r). (A.5) 

Taking the divergence of each side, we have 

(V ■m)S 3 (r) = V 2 c/ )m (r), (A.6) 

where V • {{aqH) = has been used. Comparing it with 

V 2 (V •m)Go(r) = -(V • m)<5 3 (r), (A.7) 

which is obtained from V 2 G*o(r) = — S 3 (r), we find the 
solution of (|A.6j) to be 

m (r) = -(V-m)G (r). (A.8) 

With Eq. (|A~5)) and VG = -Gi we have Eq. ©: 

H(r) = -(m- V)Gi(r) +m6 3 (r), (A.9) 

which can be confirmed to satisfy the Maxwell equations. 
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